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Abstract: Mining activity affects the vegetation and soils of the ecosystems. However, the effects of mining 
activity on saxicolous lichen communities are less concerned. Thus, the aim of this study was to characterize 
saxicolous lichen communities in three basins (Vis-Vis River basin, Poteros River basin, and Capillitas River 
basin) surrounding metalliferous mining projects of different types of operation and at different stages of 
exploitation. A large-scale mine (Bajo de la Alumbrera) with more than 25 a of open-pit mining located in 
the Vis-Vis River basin (CRV). A pre-exploitation mine (Agua Rica) located in the Poteros River basin (CRP), 
and a small-scale mine (Minas Capillitas) with more than 160 a of underground mining located in the 
Capillitas River basin (CAC). In each basin, species richness, cover, and frequency of lichen communities 
were measured on 40 rock outcrops. Also, explanatory variables were recorded, i.e., altitude, slope, aspect, 
vegetation cover, rock, and soil cover around the rocky area sampled. Richness and total cover of lichen 
communities were analysed using linear models, and species composition was explored using multivariate 
ordination analysis. Results showed that a total of 118 lichen species were identified. The species richness 
differed among basins and the lichen composition present in areas close to mining sites responded mainly to 
basins, altitude, and microsite variables. The lichen cover showed no difference among basins, but it 
changed under different rock and vegetation cover. It was not possible to quantify the effects of mining 
activity on species richness and composition. However, the low richness values found in the downstream 
of Minera Alumbrera could be associated with the negative impact of open-pit mining. Moreover, the 
effects of large-scale mining activity on lichen communities needs more investigation. 
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1 Introduction 


Mining is the main productive activity in the Catamarca Province, Argentina. Due to the current 
level of development and the geological mining potential of the region, this activity is of great 
economic importance for both the province and the country. Given the environmental threat it 
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poses (Moretton et al., 1996; Boamponsem et al., 2010; Panagos et al., 2013; Touceda-Gonzalez 
et al., 2017), mining activity must be continuously monitored to evaluate possible impacts on the 
air, soil, and water. 

Metalliferous mining is concentrated in the west of the Catamarca Province, where the Farallon 
Negro Volcanic Complex (FNVC) is located. Bajo de la Alumbrera is one of the world's largest 
open-pit copper and gold deposits. According to Alvarez (2002), the start-up of this project 
implied an investment of USD 1200x10 to produce an average of 6x10° t of concentrate per year 
through crushing, grinding, and flotation processes; and the concentrate contained approximately 
18x10^ t of copper and 20 t of gold. 

The FNVC is an eroded stratovolcano (Llambías, 1972). Within its caldera, there are vetiform 
and porphyry-type epithermal copper deposits. Among them, Farallón Negro and Altos de la 
Blenda have veins of quartz and carbonates containing the minerals aurum (Au), argentum (Ag), 
and manganese (Mn). Other deposits such as Agua Rica and Bajo de la Alumbrera are of porphyry 
type with associated metallic mineralization of copper (Cu), molybdenum (Mo), and Au. At 
Minas Capillitas, a vetiform type deposit, mineralization is mainly Cu, plumbum (Pb), zinc (Zn), 
with arsenic (As), antimony (Sb), Au, and Ag as accessory elements and wolfram (W), tin (Sn), 
bismuth (Bi) and germanium (Ge) as trace elements in gangue of quartz and alunite 
(Márquez-Zavalía, 2002; Márquez-Zavalía et al., 2014). 

Currently, the Farallón Negro and Alto de la Blenda projects are in exploitation phase, Bajo de 
la Alumbrera is in mine closure phase, while Agua Rica is in exploration phase. In Minas 
Capillitas, the mineralization zone was initially exploited for gold extraction, later for Cu at a 
smaller scale, and since the late 19" century and up to the present day, rhodochrosite has been the 
only mineral extracted for commercialization. 

The mining activity generates an effect on the biota that surrounds the work area, mainly on the 
flora. Within this stress, it has been found that mining activity generates changes in the phenology 
of nearby plants. Sun et al. (2022) found that this change was due to the decrease in the 
groundwater level caused by deep excavation, water pollution, and air pollution by dust. All these 
factors alter the metabolism of the plants and generate their greatest effect in the areas closest to 
the mines and decrease as they move further away. The structure of the plant community is also 
affected by mining, and in mining areas the plant cover always decreases (Yu and Zahidi, 2023), 
within that the cover and richness of trees and shrubs decrease, while the herbaceous cover 
increases (Kuffour et al., 2020). Regarding the diversity of species, Unanaonwi and Amonum 
(2017) measured the diversity of plant species in forests around mining areas, finding that near 
the mines the diversity values are low while these values increase as they move away from the 
mines. 

A lichen is a self-sustaining ecosystem formed by the interaction of an exhabitant fungus and 
an extracellular arrangement of one or more photosynthetic partners and an indeterminate number 
of other microscopic organisms (Hawksworth and Grube, 2020). Lichens are sensitive to 
atmospheric pollutants, with susceptibility varying among species (Pescott et al., 2015; Abas, 
2021; Paoli et al., 2021). Therefore, changes in abundance and species diversity of lichen 
communities can serve as an indicator of the adverse ecological effects of air pollution. For 
example, in polluted sites, most air pollution-sensitive lichen species are absent or declining, and 
lichen diversity is low. This effect can be caused by different mechanisms. For example, sulfur 
dioxide (SO?) can have a direct inhibitory (toxic) effect on lichens, reducing global lichen 
abundance and species richness in the community (Nash and Gries, 2002). Other pollutants, such 
as nitrogen, can stimulate the colonization and growth of some nitrophilous lichen species, 
altering lichen species diversity and community composition through changes in interspecific 
competitive relationships (Filippini et al., 2020). 

On the other hand, the abundance and diversity of lichen communities vary according to abiotic 
factors that operate at different scales; therefore, the effects of pollutants must be studied together 
with these other sources of variation (Lücking and Matzer, 2001). It is well-known that lichens, as 
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most organisms, respond to factors that change with altitude such as temperature and humidity 
(Baniya et al., 2010; Vittoz et al., 2010; Bassler et al., 2016; Rodríguez et al., 2017; Cleavitt et al., 
2019; Vetaas et al., 2019). At another scale and, depending on the latitude, other important 
variables are aspect (mainly North-South) and slope, since both determine the amount of solar 
radiation that, together with altitude, impact the incident solar radiation (insolation) and 
evapotranspiration (Kidron and Termina, 2010; Rodriguez et al., 2017; Costas et al., 2021; 
Rutherford and Rebertus, 2022). Thus, in the southern hemisphere, south-facing slopes are colder 
and more humid than north-facing slopes (Kórner, 1995, 2007). Another factor involved in the 
variation of lichen communities 1s the structure of the surrounding vegetation, since it impacts not 
only the availability of substrates but also the shading on the surface of the rock or soil. These 
explanatory variables of the changes in lichen community are very important and should be 
considered at different sites along the altitudinal gradient when designing a biomonitoring study 
associated with environmental disturbances. This is especially important in mountain areas, where 
environmental variables change with altitude, causing variation in richness and composition of 
species (Kórner, 2007). 

Although lichens from metal-enriched habitats have attracted much attention at the global scale 
(Osyczka and Rola, 2019; Osyczka et al., 2021; Neitlich et al., 2022), few efforts have been made 
to document lichen communities in mining areas (Bielczyk et al., 2009; Rajakaruna et al., 2011). 
In Argentina, lichen species from near mines areas have not been documented in detail. 

Studies of lichen communities in Argentina are scarce. Estrabou et al. (2010) studied areas of 
the province of Catamarca with the aim of establishing a baseline that would allow detection of 
environmental changes. The aim of the present study was to characterize the lichen-associated 
funga in areas surrounding mining projects and to interpret how species composition responds to 
changes in the natural environment (slope, aspect, altitude, rock, soil, and vegetation cover) and 
the influence of mining activities. In turn, taking into account the lack of works on lichen 
communities in the region, and particularly in areas with mining activity (Abas, 2021), we sought 
to contribute to the knowledge of the diversity of these organisms at high altitudes close to mining 
areas in northwestern Argentina. 


2 Materials and methods 


2.1 Study area 


The study area, located in the northwestern Sierras Pampeanas in Argentina, including sections of 
the FNVC and the western flank of the Sierra del Aconquija (Ramos, 1999). 

Due to the geographical factors that determine the inaccessibility of a large part of the study 
area, we conducted sampling in three areas located along the downstream of basins nearby 
metalliferous mining. These basins are associated with different types and stages of mining 
activity. Vis-Vis River basin (CRV) is close to Bajo de la Alumbrera that is a large-scale mine 
with more than 25 a of open-pit mining, and is currently in the mine closure stage. Poteros River 
basin (CRP) is close to Agua Rica that is planned to be a mine of the same scale as Bajo de la 
Alumbrera but still in pre-exploitation stage, and therefore we could consider CRP as a control 
area with no mining. Capillitas River basin (CAC) is close to Minas Capillitas that is a 
small-scale mine with more than 160 a of underground mining. 

The climate in the study area is temperate and markedly continental. According to the 
Kóepenand Geiger classification, it corresponds to the BWk type, i.e., cold desert. It is 
characterized by its aridity, thermal amplitude, and strong insolation, with precipitation 
concentrated in summer and a marked water deficit. Average temperatures in the study area is 
25.6?C in summer and 13.4?C in winter. Winds are dry and atmospheric humidity is very low 
(Paoli, 2003). There is no humidity data for the entire area, however, the weather station at the 
Alumbrera tailings dam for the period June, August, and September 2017 measured 39% average 
relative humidity (Minera Alumbrera YMAD-UTE, 2017). South-easterly, southerly, and westerly 
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winds predominate (Gonzalez-Bonorino, 1972). Meteorological data were recorded in Andalgala 
and Belén, the two most important localities, close to the study area. The average precipitation 
recorded at the Andalgala weather station (70 a of data) is 310 mm/a, with peaks in 1923 and 
1977, and droughts in 1941 and 1950. In the Andalgala River basin, near to Agua Rica, average 
annual precipitation is 560 mm, with a maximum of 920 mm and a minimum of 405 mm (Comba, 
2017). In the town of Belén, average annual precipitation is 244 mm and average temperature is 
25.6°C in summer and 13.4°C in winter (Gonzalez-Bonorino, 1972; Paoli, 2003). 

This area includes several phytogeographic provinces: Monte, Prepuna, and Chaco Serrano. 
According to the species present and the physiognomy of the dominant vegetation, we found that 
the area surrounding Minera Alumbrera and CRV has the botanical affinities to be included in the 
Monte phytogeographic Province (Morlans, 1995; Karlin et al., 2017). The Prepuna mainly 
occupies the slopes and high foothills of the Sierras Pampeanas, such as the southeast of the 
Sierra de Aconquija (Andalgala department) and most of the Sierra de Capillitas (Andalgala, 
Belén, and Santa Maria). CAC is located in this sector. The slopes and peaks of the mountain 
ranges between 600 and 3000 m a.s.l., belong to the Chaco Serrano. In Andalgalá, the Chaco 
Serrano includes sectors of the mountain systems of the Sierra de Capillitas and the western 
slopes of the Sierra de Aconquija, where the Agua Rica mining project and CRP are located 
(Morlans, 1995). 


2.2 Site selection 


The study was conducted between 2015 and 2017, at 8 sampling sites selected in each of the 
studied basins. Both CRV and CRP basins present fairly wide altitudinal gradients. CRV sites are 
located between 1340 and 2070 m a.s.l., and CRP sites between 1850 and 2300 m a.s.l., with the 
lowest altitudes being those furthest from the centre of mining operations (Fig. 1). On the other 
hand, CAC sites were located at the same altitude, and were distributed to the east and north of 
the mine (Fig. 1). 

2.3 Survey of lichen frequency and coverage 


At each sampling site, 5 rock outcrops (face with an aspect between 160 and 200 degrees and a 
slope between 60 and 90 degrees) were randomly selected. These conditions correspond to the 
highest values of lichen diversity in mountain environments in the Southern Hemisphere due to 
low levels of insolation (Rodríguez et al., 2017; Costas et al., 2021). A 20 cmx20 cm grid was 
placed on each rock, the lichen species were identified and the percentages (%) of relative 
frequency and relative cover of each species were estimated (Scheidegger et al., 2002; Rodríguez 
et al., 2017). In addition, as environmental or explanatory variables, the altitude of the sampling 
sites and the following microsite variables were measured: slope and aspect of the grid and the 
percentages of vegetation, rock, and soil cover around the rocky area sampled (4 m°). 


2.4 Species identification 


Morpho-anatomical and chemical analysis was performed using routine techniques. Lichen 
substances with chemotaxonomic value were identified using recrystallization and thin-layer 
chromatography techniques (Orange et al., 2001) and following keys and descriptions for genera 
and species (Adler, 1992; Scutari, 1992; Estrabou, 1999; Estrabou et al., 2006; Filippini et al., 
2014). The nomenclature of genera and species follows Calvelo and Liberatore (2002), Lücking et 
al. (2017), and Index Fungorum (http://www.indexfungorum.org/names/names.asp). Samples with 
very small or no reproductive structures were identified to genus level or with artificial names. 
One specimen of each identified species was deposited in the LUTI herbarium and in the 
Catamarca Research and Transfer Center. 


2.5 Data analysis 


Matrices of relative frequency and cover of species and environmental variables per sampling unit 
(grid) were prepared with the data obtained. From the relative frequency matrix, richness (number 
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Fig. 1 Location of sampling sites in the study area. (a), overall of the study area including the three river basins; 
(b), Vis-Vis River basin (CRV); (c), Poteros River basin (CRP); (d), Capillitas River basin (CAC). 
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of species per rock) was calculated, while from the matrix of species cover, the total cover of 
lichen was calculated per rock. The two variables were considered univariate response variables. 

To analyze richness, we applied generalised linear models with Poisson distribution, 
considering each grid as a sampling unit. In the model, the basin was considered as the main 
factor and altitude, a covariable together with percentage of rock, vegetation and soil cover, slope, 
and aspect (sine and cosine were calculated to obtain the North-South (N-S) and the East-West 
(E-W) components, respectively). The sites were considered as a random factor. Akaike 
information criterion (AIC) was used to select the model with the best fit (Di Rienzo et al., 2017). 
The distance to the mine and the altitude could not be separated as independent variables since the 
sites at the highest altitude are the closest to the mines and the sites at the lowest altitude are the 
furthest away. The total cover of lichens was analysed using a general and mixed model with the 
same explanatory variables as those used for richness. The similarity (shared species) between 
basins was analysed using a Venn diagram. 

We analyzed the variation in the composition of saxicolous lichen species between basins and 
between sites in each basin using a canonical correspondence analysis (CCA) with the aim to 
determine grouping or separation of species according to environmental or explanatory variables, 
using the species frequency matrix and the matrix of environmental variables (altitude, N-S and 
E-W component, slope, rock, and vegetation and soil cover) by 120 grids. To obtain a clear 
ordering in this analysis, we removed species with fewer than three occurrences in the sampling 
units (McCune et al., 2002). The environmental variables that had a correlation r>0.20 with the 
principal axes were plotted, which allowed us to discard those that did not influence the first two 
axes of the multivariate analysis. 

To determine the indicator species in each basin, we carried out an indicator species analysis 
(ISA) based on the method of Dufréne and Legendre (1997). This method is very useful to detect 
and describe groups of species, indicating their affinity to particular sets of environmental 
conditions. ISA provides an indicator value (IV) for each species in each group. The analysis 
combines information on the abundance of species with the occurrence within particular groups, 
with the best indicator being the species that is always present and is exclusive to the group 
(IV=100). In turn, the Monte Carlo test was applied to estimate the statistical significance of each 
IV with a probability value of 0.05% (McCune et al., 2002). All the analyses were performed 
using the Infostat software (Di Rienzo et al., 2020) and software PC-ORD Multivariate Analysis 
of Ecological Data (McCune and Mefford, 1998). 


3 Results 


A total of 118 taxa belonging to 43 genera and 18 families were identified. Of these, 83 taxa were 
identified, 18 to the genus level and 17 to the species level. According to the growth form, 55 
crustose, 42 foliose, 9 fructiculose, 4 umbilicated foliose, 5 squamulose, 1 subfructiculose, 1 
placoid, and 1 foliose-squamulose taxa were surveyed. The most diverse genera found in the 
study area were Xanthoparmelia (14 species), Caloplaca (8), Acarospora (6), Punctelia (5), and 
Rinodina (5). Sixty-two species were found at fewer than 3 sites (Table S1). 

The basin that presented the highest number of species was CRP, with 71 species, followed by 
CAC and CRV with 62 and 41 species, respectively. Thirty-six species were present only in CRP, 
25 species only in CAC, and 15 species only in CRV. CRP and CAC shared 29 species, whereas 
CRP shared 18 species with CRV. CAC and CRV shared 21 species. Only 13 shared species were 
found in the three basins (Fig. 2). CRP showed the highest proportion of foliose lichens (47.994), 
whereas the crustose morphotype was dominant in the other two basins (51.2% in CRV and 
51.2% in CAC, respectively) (Table S1). 

Richness differed significantly among basins, with CAC showing the highest average values 
(11.05), followed by CRP (7.08) and CRV (4.80) (P=0.0001; Fig. 3). None of the explanatory 
variables included in the model was significant for species richness (P>0.05), except for altitude, 
which was found to have an influence on the significantly higher values of CAC (P=0.0211). In 
the model that had the best fit (AIC), basin was the main factor, altitude was a covariate nested in 
the basin, and site was a random factor (Table 1). 
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Fig.2  Venn diagram showing the number of lichen species shared among the three basins close to metalliferous 
mines. CRV, Vis-Vis River basin; CRP, Poteros River basin; CAC, Capillitas River basin. 


Table 1 Best fit models for species richness and total lichen cover 


Index Variable Deviance P value 
Richness Basin 64.13 0.0001 
Basin>Altitude 9.72 0.0211 

F value 
Total lichen cover Rock cover 5.39 0.0226 
Vegetation cover 4.14 0.0448 


Note: Deviance corresponds to the generalised linear model with Poisson distribution applied to richness, and F value corresponds to the 
general linear model applied to the total lichen cover. 


Total lichen cover did not show a significant difference among basins (P>0.05), being almost 
the same in CAC (50.84 (+14.04)) and CRV (51.12 (£25.59)), and slightly lower in CRP (48.27 
(+21.20)) (Fig. 4). However, some microsite variables influenced total lichen cover in each site, 
being higher in sites with higher rock cover (P=0.0226) and lower with higher vegetation cover 
(P=0.0448) around them. 

The result of CCAs, using the relative frequency matrix including the three basins (Fig. 5), 
shows that the units are ordered by altitude and, to a lesser extent, soil and rock covers, with a 
canonical correlation coefficient between axis 1 and altitude of —0.832, between axis 2 and soil 
cover of —0.508, and between axis 3 and rock cover of 0.595. The points representing CAC show 
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Fig. 3 Richness of saxicolous lichens in the three basins associated with mining activity. CRP, Potreros River 


basin; CAC, Capillitas River basin; CRV, Vis-Vís River basin. Black boxes show mean values. Black boxes show 
mean values. Bars are standard deviations. 
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Fig. 4 Cover of saxicolous lichens in the three basins associated with mining activity. CRP, Potreros River basin; 
CAC, Capillitas River basin; CRV, Vis-Vis River basin. Black boxes show mean values. Bars are standard 
deviations. 


a tendency to separate from CRP and CRV points, the former being to the left of axis 1. These 
points are associated with species such as Umbilicaria haplocarpa Nyl., Caloplaca aff. sonorae 
Wetmore, Teloschistes nodulifer (Nyl.) Hillmann, Caloplaca aff. altoandina (Malme) Zahlbr., 
Candelariella vitellina (Hoffm.) Mill. Arg., and Lecidella aff. granulosula (Nyl.) Knoph & 
Leuckerts, also with a high IV in the ISA (Table 2). CCA had a total inertia in the species data of 
11.33 and a total variance of 8.4% explained by the three axes. 

The ISA showed that 34 species had significant IVs for the different basins (P<0.05). All the 
groups (basins) had different indicator species that were specific to those environmental 
conditions (Table 2). Eighteen species were indicators of CAC, 10 species were indicators of CRP, 
and 6 species were indicators of CRV. 


4 Discussion 


Most studies on lichens associated with mining focus on the restoration of lichen communities 
after mining exploitation (Abas, 2021). As far as we know, studies about metalliferous mining and 
the lichen community biomonitoring approach are scarce. The number of lichens surveyed in this 
work (118) is high compared with previous ecological studies carried out in central-western 
Argentina. Previous studies of saxicolous lichens in central and north-western Argentina reported 
107 species in the mountains of Cordoba (Rodriguez et al., 2017) and 58 species in the mountains 
of La Rioja (Costas et al., 2021). 

Results of the species richness and composition analyses showed great differences in the lichen 
communities among the three basins, suggesting environmental differences. The role of altitude in 
lichen composition and abundance (Rodriguez et al., 2017) explains the highest values of species 
richness per sampling unit found in CAC, since this basin has a more temperate climate than the 
other two basins, which favours the development of lichen communities. Similarly, Costas et al. 
(2021) recorded the highest richness values at 2897 m a.s.l., coinciding with the altitudinal belt of 
the sampling sites in CAC. However, although CRP has lower lichen average richness values per 
sampling unit than CAC, it has the highest number of species (71; Fig. 1), probably due to the 
wide altitudinal gradient between sampling sites. Thus, an increase in species richness was 
observed with increasing altitude, suggesting that changes in environmental conditions with 
increasing altitude favour the development of species that are not present at lower altitudes. 

CRV is the area with the lowest richness and diversity. Taking into account that it has an 
altitudinal gradient similar to that of CRP, altitude is not the main factor causing differences in 
lichen richness between these two basins. On the other hand, CCA (Fig. 3) showed that CRV 
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Fig. 5 Canonical correlation analysis (CCA) using the relative frequency of lichen species. Grey arrows show 
the correlation of each axis with altitude and soil cover. CRP, Potreros River basin; CAC, Capillitas River basin; 
CRV, Vis-Vis River basin; A.alto, A. altoandina; A.chry, Acarsopora chrysops (Tuck. H.Magn; A.lore, 
Acarospora lorentzii (Müll. Arg.) Hue; A.xant, Acarospora xanthopahana (Nyl.) Jatta; B.spur, Buellia spuria 
(Schaer.) Anzi; B.sulf, Buellia sulphurea Malme; C.aff.al, Caloplca aff. altoandina (Malme) Zahlbr.; C.ama3, 
C.amer 3; C.atro, Caloplaca atroflava (Turner) Mong.; C.conc, Candelaria concolor (Dicks.) Arnold; C.nubi, 
Culbersonia nubile (Moberg) Essl.; C.ochr, Caloplaca ochraceofulva (Müll. Arg.) Jatta; C.side, Caloplaca aff. 
sideritis (Tuck.) Zahlbr.; C.sono, C. sonorae; C.vite, C. vitellina; F.hays, Flavoparmelia haysomii (C.W. Dodge) 
Hale; Fla.sp, Flavoplaca sp; L.argo, Lecanora argopholis (Ach.) Ach.; L.carp, Lecidella carpathica Kórb.; L.disp, 
Lecanora dispersa (Pers.) Róhl.; L.gran, L. granulosula; L.mura, Lecanora muralis (Schreb.) Rabenh.; Lep.sp, 
Lepraria sp.; N.pulc, Nephroma pulchella (Borrer) Nyl.; P.aipo, Physcia aipolia (Ehrh. ex Humb.) Fürnr.; P.badi, 
Protoparmelia badia (Hoffm.) Hafellner; P.breu, Pyrenula aff. breutelii (Müll. Arg.) Aptroot; P.cand, 
Placomaronea candelarioides Räsänen; P.dist, Psiloparmelia distinta (Nyl.) Hale; P.micr, Punctelia microsticta 
(Müll. Arg.) Krog; P.peer, Punctelia perreticulata (Rásánen) G. Wilh. & Ladd; P.ponc, PAyscia poncinsii Hue; 
P.punc, Punctelia punctilla (Hale) Krog; P.reti, Parmotrema reticulatum (Taylor) M. Choisy; P.sema, Punctelia 
semansiana (W.L. Culb. & C.F. Culb.) Krog; P.stic, Punctelia stictica (Delise ex Duby) Krog; R.cela, Ramalina 
celastri (Spreng.) Krog & Swinsc.; R.conf, Rinodina confragosa (Ach.) Kórb.; R.disp, Rhizocarpon disporumn 
(Nàgeli ex Hepp) Müll. Arg.; R.inca, Ramalina incana Kashiw.; R.long, Rinodina longisperma Matzer & H. 
Mayrhofer; T.nodu, Teloschistes nodulifer (Nyl.) Hillmann; T.sedi, Toninia sedifolia (Scop.) Timdal; U.calv, 
Umbilicaria calvescens Nyl.; U.duri, Usnea durietzii Motyka; U.ploy, Umbilicaria polyphylla (L.) Baumg; X.cord, 
Xanthoparmelia cordillerana (Gyeln.) Hale; X.fari, Xanthoparmelia farinosa (Vain.) T.H. Nash, Elix & J. Johnst.; 
X.hass, Xanthomendoza hasseana (Räsänen) Sechting, Kàárnefelt & S.Y. Kondr.; X.loxo, Xanthoparmelia loxodes 
(Nyl.) O. Blanco, A. Crespo, Elix, D. Hawksw. & Lumbsch; X.micr, Xanthoparmelia microspora (Müll. Arg.) 
Hale; X.punc, Xanthoparmelia punctulata (Gyeln.) Hale; X.sant, Xanthoparmelia santessonii T.H. Nash & Elix. 


is separated from CRP, being influenced by the soil cover. In CRV, where soil cover is higher than 
in CRP, vegetation cover is poorer, which in turn influences lichen composition through greater 
light availability on rocks. In a work carried out in mountain systems similar to those studied here, 
as altitude increases, the proportion of crustose lichens increases (Costas et al., 2021). This 
phenomenon occurs in CAC, which is the highest basin and the one with the highest number of 
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Table 2 Indicator Species Analysis (ISA) according to basins and applied to presence-absence lichen data 


Species Basin IV GF P 
Xanthoparmelia punctulata (Gyeln.) Hale CAC 25.0 F 0.05 
Protoparmelia badia (Hoffm.) Hafellner 12:5 [o 0.01 
Ramalina incana Kashiw. 9.7 Fr 0.05 
Xanthoparmelia cordillerana (Gyeln.) Hale 24.6 F 0.00 
Lecidella aff granulosula (Nyl.) Knoph & Leuckert 29.4 [o 0.00 
Punctelia perreticulata (Rásánen) G. Wilh. & Ladd 21.2 F 0.01 
Umbilicaria haplocarpa Nyl. 59.7 FU 0.00 
Psiloparmelia distincta (Nyl.) Hale 21.9 F 0.00 
Xanthoparmelia microspora (Müll. Arg.) Hale 24.0 F 0.00 
Xanthomendoza hasseana (Räsänen) Sochting, Kárnefelt & S. Y. Kondr. 17.1 F 0.00 
Caloplaca aff sonorae Wetmore 67.5 C 0.00 
Teloschistes nodulifer (Nyl.) Hillmann 27.9 Fr 0.00 
Caloplaca aff. altoandina (Malme) Zahlbr. 52:5 C 0.00 
Candelariella vitellina (Hoffm.) Müll. Arg. 33.3 [o 0.00 
Lecidella carpathica Kórb. 28.1 C 0.00 
Xanthoparmelia loxodes (Nyl.) O. Blanco, A. Crespo, Elix, D. Hawksw. & Lumbsch 17.5 F 0.00 
Rinodina confragosa (Ach.) Kórb. 15.4 [o 0.01 
Buellia sulphurea Malme 17.6 C 0.03 
Caloplaca atroflava (Turner) Mong. CRP 30.9 C 0.00 
Buellia spuria (Schaer.) Anzi 51.9 [o 0.00 
Lecanora argopholis (Ach.) Ach. 39.6 C 0.00 
Punctelia borrerina (Nyl.) Krog 17.3 F 0.00 
Ramalina celastri (Spreng.) Krog & Swinscow 20.0 Fr 0.00 
Usnea durietzii Motyka 21.9 Fr 0.00 
Umbilicaria calvescens Nyl. 20.0 FU 0.00 
Normandina pulchrella (Borrer) Nyl. 10.0 E 0.03 
Physcia aipolia (Ehrh. ex Humb.) Fürnr. 12.5 F 0.01 
Rhizocarpon disporum (Nàgeli ex Hepp) Müll. Arg. 12.2 C 0.03 
Acarospora lorentzii (Müll. Arg.) Hue CRV 17.8 C 0.02 
Caloplaca ochraceofulva (Müll. Arg.) Jatta 27.9 C 0.00 
Lecanora muralis (Schreb.) Rabenh. 12.0 P 0.03 
Punctelia punctilla (Hale) Krog 24.6 F 0.00 
Rinodina longisperma Matzer & H. Mayrhofer 29.2 C 0.00 
Caloplaca aff americana (Malme) Zahlbr. 29.3 C 0.00 


Note: CAC, Capillitas River basin; CRP, Potreros River basin; CRV, Vis-Vis River basin; IV, indicator value; GF, growth form; F, 
foliose; Fr, fruticulose; FU, foliose umbilicate; C, crustose. E, squamulose; P, placoid. 


crustose species. Since CRP and CRV have similar altitudinal gradients, the lower percentages of 
foliose species found in CRV than in CRP could be related to the intrinsic environmental 
heterogeneity of each basin (Kórner and Spehn, 2002). However, since species diversity and 
growth form (crustose, foliose, and fruticose) are affected by atmospheric pollution 
(Gunawardana et al., 2021), the low proportion of foliose species found in CRV could be related 
to the open-pit mining activity of Bajo de la Alumbrera. 

Lichen cover values in the three basins were about 50%, with no significant differences among 
basins. However, this value is influenced by microsite variables (vegetation and rock covers), as 
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observed by Rodríguez et al. (2017) for lichen communities in central Argentina. 

The analysis of indicator species showed a close relationship between altitude and the 
occurrence of those species in different basins. In CRV, A. xanthophana and P. candelarioides 
were associated with higher altitude, which is consistent with the habitat range observed in 
taxonomic descriptions for this species (Knudsen et al., 2008; Westberg et al., 2009; Knudsen and 
Flakus, 2016). On the other hand, in CRP, Xanthoparmelia farinosa (Vain.) T.H. Nash, Elix & J. 
Johnst. was associated with lower altitude, while P. perreticulata and U. durietzzi were indicators 
of high-altitude sites. This pattern coincides with the habitat range observed in central Argentina 
(Estrabou, 1999; Rodríguez et al., 2011). In this sense, the highest number of indicator species 
found in CAC shows the high specialization of this altitude (3000 m a.s.l.). The species with 
IV»50 in this area were C. sonorae, U. haplocarpa, and C. altoandina. U. haplocarpa is endemic 
to the central Andes (2500-4400 m a.s.1.) (Hestmark, 2010). 


5 Conclusions 


Saxicolous lichen communities were studied with the aim of contributing to the knowledge of the 
diversity of lichens that grow in high-altitude mining areas in the west of the Catamarca province, 
Argentina. The characteristics of the lichen communities responded to environmental factors 
(mainly altitude and weather conditions) and microsite variables. Therefore, it was not possible to 
attribute the differences in species richness and composition among basins and sites to mining 
activity. However, the very low richness values found in CRV suggest that open-pit copper mining 
could influence lichen composition. More in-depth ecological studies analysing the lichen 
communities inside and outside the mine in a narrow altitudinal range will allow us to more 
clearly determine the impact of mining activities on lichens. Finally, taking into account the lack 
of works on lichen communities in the province, and particularly those associated with mining, 
these results lay the groundwork for future in-depth studies on the impact of mining activity on 
lichen species composition in the western region of Catamarca. 
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Appendix 
TableS1 List of lichen taxa identified by basin, taxonomic family, and growth form 
Taxon = Family cane 
CRP CAC CRV form 
Acarospora aff. altoandina H. Magn. F Acarosporaceae C 
Acarospora aff. obnubila H. Magn. + Acarosporaceae C 
Acarospora boliviana H. Magn. H Acarosporaceae C 
Acarospora lorentzii (Müll. Arg.) Hue + + + Acarosporaceae C 
Acarosporoide negro + Acarosporaceae C 
Acarospora xanthophana (Nyl.) Jatta + + Acarosporaceae C 
Acarospora chrysops (Tuck.) H. Magn. + + Acarosporaceae C 
Buellia aff. umbrina Malme i Caliciaceae C 
Buellia spuria (Schaer.) Anzi + + + Caliciaceae C 
Buellia sulphurea Malme + + Caliciaceae C 
Caloplaca aff. americana (Malme) Zahlbr. + + Teloschistaceae C 
Caloplaca aff. sideritis (Tuck.) Zahlbr. * Teloschistaceae C 
Caloplaca aff. sonorae Wetmore + Teloschistaceae C 
Caloplaca aff. subsquamosa (Müll. Arg.) Zahlbr. + Teloschistaceae C 
Caloplaca amarillo + Teloschistaceae C 
Caloplaca atroflava (Turner) Mong. + + Teloschistaceae C 
Caloplaca ochraceofulva (Müll. Arg.) Jatta + + + Teloschistaceae C 
Caloplca aff. altoandina (Malme) Zahlbr. * Teloschistaceae C 
Candelaria concolor (Dicks.) Arnold + + Lecanoraceae F 
Candelaria fibrosa (Fr.) Müll. Arg. + Lecanoraceae FẸ 
Candelariella vitellina (Hoffm.) Müll. Arg. + + Candelariaceae C 
Cladonia pyxidata (L.) Hoffm. + Cladoniaceae F 
Crustoso amarillo 1 + - C 
Crustoso amarillo 2 + - C 
Crustoso amarillo 3 * - C 
Crustoso blanco 1 + - C 
Crustoso blanco 2 + - C 
Crustoso blanco 3 + - C 
Crustoso blanco con peritecios + - C 
Crustoso gris verdoso * - C 
Crustoso negro F - C 
Crustoso sp. 1 + = C 
Crustoso verde - C 
Culbersonia nubila (Moberg) Essl. + Physciaceae F 
Diploschistes euganeus (A. Massal.) J. Steiner F Thelotremataceae C 
Diploschistes bisporus (Bagl.) J. Steiner + Thelotremataceae C 
Diploschistes scruposus (Schreb.) Norman * Thelotremataceae C 
Diplotomma hedinii (H. Magn.) P. Clerc & Cl. Roux&nbsp + Caliciaceae C 
Escuamuloso marrón + - E 


To be continued 
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Continued 
Taxon os Family Groni 

CRP CAC CRV form 
Flavoparmelia haysomii (C.W. Dodge) Hale + + + Parmeliaceae F 
Flavoparmelia papilosaun (Lynge ex Gyeln.) Hale + Parmeliaceae E 
Flavoparmelia soredians (Nyl.) Hale m Parmeliaceae F 
* S (Vondrák, Říha, Arup & Sechting) " Teloschistaceae C 
Flavoplaca sp. a3 Teloschistaceae C 
Heterodermia albicans (Pers.) Swinscow & Krog WE Physciaceae F 
Heterodermia diademata (Taylor) D.D. Awasthi + Physciaceae F 
Heterodermia lutescens (Kurok.) Follmann + Physciaceae SFr 
Heterodermia obscurata (Nyl.) Trevis. + Physciaceae F 
Hyperphyscia syncolla (Tuck. ex Nyl.) Kalb + Physciaceae F 
Lecanora argopholis (Ach.) Ach. + + Lecanoraceae C 
Lecanora dispersa (Pers.) Róhl. + Lecanoraceae C 
Lecidella aff. granulosula (Nyl.) Knoph & Leuckert + + + Lecanoraceae C 
Lecidella carpathica Kórb. * * Lecanoraceae C 
Leprarioide sp. * Stereocaulaceae C 
Leptogium austroamericanum (Malme) C.W. Dodge + Collemataceae F 
Leptogium hypotrachynum Müll. Arg. + Collemataceae F 
Leptogium phyllocarpum (Pers.) Mont. + Collemataceae F 
Microfolioso marron + - F 
Normandinaephroma pulchella (Borrer) Nyl. + Nephromataceae E 
Paraparmelia sp. + Parmeliaceae F 
Parmotrema consors (Nyl.) Krog & Swinscow F Parmeliaceae F 
Parmotrema muelleri (Vain.) O. Blanco, A. Crespo, Divakar, ; 
Elix & Lumbsch + Parmeliaceae F 
Parmotrema reticulatum (Taylor) M. Choisy + + Parmeliaceae F 
Phaeophyscia hirsuta (Mereschk.) Essl. + Physciaceae F 
Physcia aipolia (Ehrh. ex Humb.) Fürnr. + Physciaceae F 
Physcia poncinsii Hue + + Physciaceae F 
Physcia tribacia (Ach.) Nyl. + Physciaceae F 
Placidium ruiz-lealii (Räsänen) Breussruiz-lealli * Verrucariaceae E 
Placomaronea candelarioides Räsänen + + + Lecanoraceae FU 
Protoparmelia badia (Hoffm.) Hafellner + Parmeliaceae C 
Psiloparmelia distincta (Nyl.) Hale + + Parmeliaceae F 
Punctelia borrerina (Nyl.) Krog + + Parmeliaceae F 
Punctelia perreticulata (Räsänen) G. Wilh. & Ladd + + + Parmeliaceae F 
Punctelia punctilla (Hale) Krog F $ + Parmeliaceae F 
Punctelia semansiana (W.L. Culb. & C.F. Culb.) Krog + + Parmeliaceae F 
Punctelia stictica (Delise ex Duby) Krog + + + Parmeliaceae F 
Pyrenula aff. breutelii (Müll. Arg.) Aptroot + + Pyrenulaceae C 
Ramalina celastri (Spreng.) Krog & Swinsc. F Ramalinaceae Fr 
Ramalina incana Kashiw. + + Ramalinaceae Fr 


To be continued 
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Continued 
Taxon os Family owe 
CRP CAC form 
Rhizocarpon disporumn (Nàgeli ex Hepp) Müll. Arg. + + Rhizocarpaceae C 
Rinodina confragosa (Ach.) Körb. * * Physciaceae C 
Rinodina longisperma Matzer & H. Mayrhofer * zd Physciaceae C 
Rinodina oxydata (A. Massal.) A. Massal. + + Physciaceae C 
Rinodina peloleuca (Nyl.) Müll. Arg. + Physciaceae C 
Rinodina sp. * Physciaceae C 
Staurothele areolata (Ach.) Lettau + + Verrucariaceae C 
Teloschistes chrysophthalmus (L.) Th. Fr. + Teloschistaceae Fr 
Teloschistes nodulifer (Nyl.) Hillmann + Teloschistaceae Fr 
Toninia sedifolia (Scop.) Timdal Ramalinaceae 
Toninia aff submexicana B. de Lesd. Ramalinaceae 
ee smaragdula (Wahlenb. ex Ach.) Cl. Roux Acarosporaceae C 
Umbilicaria calvescens Nyl. F Umbilicariaceae FU 
Umbilicaria haplocarpa Nyl. + + Umbilicariaceae FU 
Umbilicaria polyphylla (L.) Baumg. + + Umbilicariaceae FU 
Usnea amblyoclada (Müll. Arg.) Zahlbr. Parmeliaceae Fr 
Usnea cirrosa Motyka + Parmeliaceae Fr 
Usnea dasaea Stirt. + Parmeliaceae Fr 
Usnea durietzii Motyka + + Parmeliaceae Fr 
Usnea sp. F Parmeliaceae Fr 
Verrucarioide sp. + Verrucariaceae C 
e a hasseana (Räsänen) Søchting, Kärnefelt & " E Parmeliaceue F 
Xanthoparmelia cordillerana (Gyeln.) Hale F Parmeliaceae F 
Xanthoparmelia farinosa (Vain.) T.H. Nash, Elix & J. Johnst. + + Parmeliaceae F 
Xanthoparmelia ferraroiana T.H. Nash, Elix & J. Johnst. F Parmeliaceae F 
Xanthoparmelia fumarprotocetrarica (Elix & J. Johnst.) Elix + Parmeliaceae F 
il ua ese n a (Nyl.) O. Blanco, A. Crespo, Elix, D. + Parmeliaceue F 
Xanthoparmelia mahuiana T.H. Nash & Elix Parmeliaceae F 
Xanthoparmelia microspora (Müll. Arg.) Hale + + Parmeliaceae F 
Xanthoparmelia mougeotii (Schaer. ex D. Dietr.) Hale + Parmeliaceae F 
Xanthoparmelia punctulata (Gyeln.) Hale + + Parmeliaceae F 
Xanthoparmelia santessonii T.H. Nash & Elix H Parmeliaceae F 
Xanthoparmelia subulcerosa T.H. Nash & Elix + Parmeliaceae F 
Xanthoparmelia ulcerosa (Zahlbr.) Hale + Parmeliaceae E 
Xanthoparmelia villamilianus T.H. Nash & J. Johnst. + Parmeliaceae F 


Note: +, presence; C, crustose; F, foliose; Fr, fruticose; SFr, sub-fruticose; E, squamulose; FU, foliose umbilicate; P, placoid. 


